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The influence of nonperturbative flelds on instantons in quantum chromodynamics is studied. 
Nonperturbative vacuum is described in terms of nonlocal gauge invariant vacuum averages of 
gluon field strength. Effective action for instanton is derived in bilocal approximation and it is 
demonstrated that stochastic background gluon fields are responsible for infra-red (IR) stabilization 
of instantons. Dependence of characteristic instanton size on gluon condensate and correlation 
length in nonperturbative vacuum is found. It is shown that instanton size in QCD is of order of 
0.25 fm. Comparison of obtained instanton size distribution with lattice data is made. 



I. INTRODUCTION 



Instanton is the first explicit example of nonperturbative (NP) quantum fiuctuation of gluon field in QCD. It was 
introduced in 1975 by Polyakov and coauthors [Q. An important development originated with 't Hooft's classic 
paper |^], in which he calculated the semi-classical tunneling rate. Instanton gas as a model of QCD vacuum was 
proposed in pioneer works by Callan, Dashen and Gross These topologically nontrivial field configurations 

are essential for the solution of some problems of quantum chromodynamics. Instantons allow to explain anomalous 
breaking of U{1)a symmetry and the 77' mass Spontaneous chiral symmetry breaking (SCSB) can be explained 

with the help of instanton and anti- instanton field configurations in QCD vacuum An important role of instan- 
tons in scalar and pseudoscalar channels was pointed out in paper ||^. Taking into account instantons is of crucial 
importance for many phenomena of QCD (see ^ and references therein). 

At the same time, there is a number of serious problems in instanton physics. The first is the divergence of integrals 
over instanton size p at big p. This makes it impossible to calculate instantons' contribution to some physical 
quantities, such as vacuum gluon condensate. Second, "area law" for Wilson loop can not be explained in instanton 
gas model, hence quasiclassical instanton anti-instanton vacuum lacks confinement which is responsible for hadron 
spectra. 

There were many theoretical works aimed to solve the problem of instanton instability in it's size p. To some extent 
they all were based on an attempt to stabilize instanton ensemble due to effects of interaction between pseudoparticles. 
The most popular is the model of " instanton liquid" , which was phenomenologically formulated by Shuryak ||l^ . It 
states that average distance between pseudoparticles is ^ ~ 1 fm and their average size is p ~ 1/3 fm. Thus, 
p/R ~ 1/3 and vacuum consists of well separated, and therefore not very much deformed, instantons and anti- 
instantons. Quantitatively close results were obtained by Diakonov and Petrov [pd]| . In their approach stabilization 
is related to discovered by them classical repulsion between pseudoparticles. However, further development |lj] 
revealed that instanton ensemble can not be stabilized due to purely classical interaction. Thus, the mechanism for 
the suppression of large-size instantons in the ensemble of topologically non-trivial fields is still not understood. 

On the other hand, QCD vacuum contains not only quasiclassical instantons, but also other nonperturbative fields. 
The investigation of instantons' behavior in stochastic QCD vacuum, which is parametrized by "vacuum correlators" 
(i.e. nonlocal gauge invariant averages of gluon field strength) was started in Refs. |l|, |l|l. It was shown, that in 
nonperturbative vacuum standard perturbation theory for instanton changes , which leads to "freezing" of effective 
coupling constant. It was then found in |l4[ that nonlocal interaction between large-size instanton (p > 1 fm) and 
nonperturbative stochastic background field does not lead to it's infrared inflation. 
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The aim of this work is to demonstrate that instanton can be stabihzed in nonperturbative vacuum and exist as a 
stable topologicaUy nontrivial field configuration against the background of stochastic nonperturbative fields, which 
are responsible for confinement, and to find quantitatively it's size. We evaluate nonlocal effective action for instanton 
in NP vacuum in gauge invariant way, describing vacuum in terms of "vacuum correlators", and show that in bilocal 
approximation effective action has a minimum at the definite instanton size pc ■ The value of pc is in functional 
dependence on properties of bilocal correlator {G{x)^{x,y)G{y)), i.e. it depends on two parameters: (G^) - gluon 
condensate and it's "measure of inhomogeneity " Tg - correlation length in the condensate. 



II. GENERAL FORMALISM 



The influence of NP fluctuations on instanton can be separated into two parts. First, perturbation theory undergoes 
a change in NP background, which results in the change of standard one-loop renormalization of instanton action. 
Second, a direct nonlocal interaction of instanton with NP background field exists. 

Standard euclidian action of gluodynamics has the form 

S[A] - ^1 d^xtr{F'^M]) d'xF^^M]F^^Ml (1) 

where .^^^.[A] = S^A^ — d^A^ — i[A^,A,j] is the strength of gluon field and we use the Hermitian matrix form for 
gauge fields A^{x) = gaAl{x)t'' /2 and tvf^t^ = 5"V2- We decompose A^, as 

A^^4-' + S,+a^, (2) 

where A"'^' is an instanton-like field configuration with a unit topological charge Q^ri^'"^'] — 1; Q/j is quantum 
field (expansion in reduces to perturbation theory, which in gluodynamics leads to asymptotic freedom); 
is nonperturbative background field (with zero topological charge), which can be parametrized by gauge invariant 
nonlocal vacuum averages of gluon field strength^. 

In general case effective action for instanton in NP vacuum takes the form 

where (...) implies averaging over background field B^, 



d{B)) = j dp{B)d{B) (4) 



and dp{B) is the measure of integration over NP fields, explicit form of which is not important for the following 
consideration. 

Expanding S[A\ up to the second power in one obtains Z = {Zi{B)Z2{B)) , where 

Zi(B) =e-^[^'""l y"[l?aj det(V2)exp ^^ir j d^x ^,a,f + 2iF^,[a^,a,]] (5) 

Z2{B) ^ exp{-5[A'"^' + B] + S-iA'"*^*]} (6) 

Here we use the notation A = A"^^^ + B, F^^ = F^^[A\, and V^Oj, = d^a^ — i[Afj_,a^] is a covariant derivative. 
Integration over a and B in (|^),(^) corresponds to averaging over fields that are responsible for the physics at different 
scales. Integration over takes into account perturbative gluons and describes phenomena at small distances. 
Averaging over (formally interaction with gluon condensate) accounts for phenomena at scales of confinement 
radius. Therefore it is physically clear that averaging factorizes and one obtains 

Z^{Z,iB)){Z2iB)) (7) 



In operator product expansion method and in QCD sum rules nonperturbative field is characterized by a set of local gluon condensates 
(G2>, (G3), ... 



It should be noted that in the hmit of infinite number of colors iVc — > oo factorization is exact Z{Nc oo) = 
{Zi{B)) {Z2{B)) . This allows us to write effective action of instanton in NP vacuum as a sum of two terms, "pertur- 
bative" and " nonperturbative" : 

5eff [A-''*] = [A-''*] + 5 J [^"-'] (8) 

^P^[A'-*] = -ln(Zi(S)) (9) 
5,NPj^i„.t] ^ _ ln(Z2(B)) - - In (exp{-5[A'""' + B] + (10) 

III. INTERACTION BETWEEN INSTANTON AND NONPERTURBATIVE VACUUM GLUON FIELDS 

The general expression for onc-instanton field configuration has the well known form 

- ^^'^''t^^^^f (^^^) ' (11) 

where matrix R!''^ ensures embedding of instanton into SU{Nc) group, h — 1,2, . . .iV^ — 1; /3 = 1,2,3, 77"^ are 
't Hooft symbols. R!"^ satisfied relations pbc j^tp j^cj ^ ePi&R^s^ R''I^R''f ^ 5'^'' and pbc j^aa j^bf3 j^c, ^ ^ap^ ^ 
singular gauge profile function f{z) satisfies boundary conditions /(O) = 1, /(oo) = and the classical solution has 
the form 



The probability to find an instanton is determined by the classical action functional 5*01 [A] for the solution (pT|) , ([L2[) ; 
that is 

w ^ exp{-5ci[Af *]} - exp{-87rVgo'} (13) 

The preexponential factor was calculated in The result for one- instanton contribution to the QCD partition 
function is = J dn{p, xq, R), where dn is the differential instanton density 

dn{p,XQ,R) [dR\d'^xo^do{p) (14) 
4.6exp{-1.687Vc} / Stt^ V'^" ( Stt"^ 



In the two- loop approximation in gluodynamics, the coupling constant g^{p) is given by 

blnf^) +^lnlnf^) +0(l/ln(l/pA)), b = ^N,, = ^N, (16) 



S(2(p) ypAj b \pA 

A = Apv corresponds to the Pauli-Villars regularization scheme, A ~ 200 MeV. 

Callan, Dashen and Gross ^] were first to show that in a constant gauge field instanton behaves as a colored four- 
dimensional dipole. Shifman, Vainshtein and Zakharov [|l^ generalized this result to the case of interaction between a 
small-size instanton, p < 0.2 fm, and nonperturbative long-wave fluctuations that are described by the local vacuum 
condensate (G^). As a result, do{p) is replaced by dcs{p), where 

dM « (Ap)" (1 + J^Y^^P' + •••) (17) 

Thus, we arrive at the well-known problem of the infrared inflation of instantons in p. 

It was shown in |l3| [l^ that in NP vacuum standard perturbation theory for instantons changes, which results 
in "freezing" of effective coupling constant. The perturbative part of effective instanton action in stochastic vacuum 
S'^fj[A'"'^*] was shown to be 



^eyp)-^nl/^^ (18) 



Here ~ 0.75mo++ ^ IGeV, where 0"^"*" is the hghtest glueball. It follows from (^8|) that for small-size instantons 
a standard perturbative result S^g{p ^ l/m*) = 87r^/g^(p) is recovered, and for large-size instantons S^g{p > 
1/m*) — > const. Correspondingly, factor ^ (Ap)'' in ( p^ ) ceases to grow with the increase of p. Next, large-size 
instanton p > po = (192/(G^})^/'* ~ 1 fm was considered in It's field is weak compared to characteristic field 
strengths in giuon condensate (G^), and it was considered as a perturbation in [l4[ |. It was demonstrated that due to 
instanton's nonlocal interaction with giuon condensate there is no infrared inflation for p > pQ. 

We consider effect of NP fields on instanton not assuming that it's field is weak (for arbitrary p), i.e. we evaluate 
(Z2). In this work we make use of the method of vacuum correlators, introduced in works of Dosch and Simonov [^ . 
NP vacuum of QCD is described in terms of gauge invariant vacuum averages of giuon fields (correlators) 

A^i^i...p„i.„ = (trG^i^,(xi)$(a;i,a;2)Gp2i'2(x2)...G^„i.„(a;„)<I>(x„,a;i)), 

f - \ . 

where G^j/ is giuon field strength, and ^{x,y) — Pexp i ^ Bf^^dz^ is a parallel transporter, which ensures gauge 

^ ^ /. 

invariance. In many cases bilocal approximation appears to be sufficient for qualitative and quantitative description 
of various physical phenomena in QCD. Moreover, there are indications that corrections due to higher correlators are 
small and amount to several percent . 

The most general form of bilocal correlator, which follows from antisymmetry in Lorentz indices, is given by 



—^{^fipSi,a - SpaS^p) H ^{Uf^npS^cr + ni,n„Sp,p - npU^S^p - n,ynpSpa)^ , (19) 

where G^^{x,xo) = ^{xq, x)Gp^{x)^{x, xq), = z^/\z\ = {x - y)p/\x - y\ is the unit vector, (G^) = (ff^G°j,G°^) 
and, as it follows from normalization, D{0) + D{Q) — 1. 

Functions D{z) and D{z) have both perturbative and nonperturbative contributions. We will consider only nonper- 
turbative part, because perturbative one has already been taken into account in S^g. Most data about nonperturbative 
components of these functions come from numerical simulations on lattice. Giuon condensate (G^) is also determined 
from lattice calculations, but there exists a widely used estimate for this value based on charmonium spectrum analysis 
and QCD sum rules (G^) ~ O.SGeV"' [|8|. According to the lattice data D{z) and D{z) are exponentially decreasing 
functions D{z) = Aq exp(— z/Tq), D{z) = Aiz exp(—z /Tg) iTg, where Tg is the gluonic correlation length, which was 
measured on the lattice ||l^ and estimated analytically to be Tg ^ 0.2 fm. Besides, according to lattice measure- 
ments Ai <^ A(i [Ai ^ v4o/10). Lattice data from paper of Di Giacomo |^ are presented in Table |[ SU[i) full stands 
for chromodynamics with 4 quarks, while SU (2) and SU (3) quenched mean pure SU (2) and SU (3) gluodynamics, 
respectively. Note, that (G^) — 0.87 GeV* for SU{i) full is in good agreement with more recent calculations based 
on QCD sum rules (22|. 

TABLE I: Lattice data for bilocal correlator 





(G^), GeV* 


Tg, fm 


SU{2) quenched 


13 


0.16 


5J7(3) quenched 


5.92 


0.22 


SU{i) full 


0.87 


0.34 



To evaluate S^i defined in (0) we use the cluster expansion, which is well known in statistical physics ]23|] : 

(exp(^))=exp(^^), (20) 



where {x) = {{x))-, {x^) = {{x^]) + {xf; {x^) = ((x^)) + 3{x){{x^)) + (x)^; 
Next, we modify expression ([lO|) for 5*^^ by adding S[B] to the exponent: 



-'cff 

_5.NPj^i„st] ^ _ ln(exp{-S'[A'"''* + B] + + S[B]}) (21) 



In bilocal approximation this is equivalent to a change of normalization of partition function^, which is not important 
for the following consideration. 

Nonperturbative part of the effective instanton action S^g takes the form 



S^i = (S'[A'"''* + B]- S[B] ~ S'[A'"''*]) - 

i (((S-i^'"^* + B]- S[B] - - (S-iA'"^' + B]- S[B] - fi-iA'"^*])^) + 



(22) 



By using Fock-Shwinger gauge x^Afj^ = Xf_,Bf^ — (instanton field A^^^^ satisfies it owing to the properties of 't Hooft 
symbols) and taking into account that 



+ 



2F^,[A"-*]G^,[i?] - 4i(F^,[A'-'] + G^,[S])[4-\i3,]} 



one gets 



(23) 



^cS — 'S'dia + T^iSdia + "^para + S'l +52, 



where in bilocal approximation 



1 



•S'dia = / '^''a^(tr 



{[A^,B,] - [A,,B^] 



(24) 

(25) 
(26) 
(27) 
(28) 

We use notations S'dia (diamagnetic) and S'para (paramagnetic) for contributions (Eq}and (26) into interaction of 
instanton with background field. Physical motivation for this is discussed in detail iii [|l4[ . For the sake of simplicity 
we drop further the index 'inst' for instanton field, A™^*' = A^. 

1 

In Fock-Shwinger gauge B^{x) — x^ J adaG,yfj_{ax), and that allows us to substitute vacuum averages in (|25|)-(|2^) 



with correlators. 

The general form of the instanton field is 



5para = -^ J d^xd^ [F^,{x)G ^,{x)) tv {Fp„{y)G p„{v))) 
Si = ^ j d^xd^yixY (F^,[A^,S,])^tr {FpMp,Ba\), 
S2^^ [ d^xd^y Uy {F^,G^,)^iY {Fp„[Ap,B„])^^ 



Ap{x) = <^{x,XQ)Ai''^{x - xo)<^{xo,x), 



(29) 



where $(x,y) = Pexp yi jB^^dz^j is a parallel transporter, and A^™^[x — xq) is instanton field in a fixed (sin- 
gular) gauge (|Tll),(|l2|). Inserting (29) into (^-(p8|) and using representation ([l9|) for gauge invariant condensate 



Indeed, starting from one finds that in bilocal approximation ( |24[ ) gets additional terms of the type (G^) / d'^x, which do not depend 
on ^4'°"* and therefore result only in renormalization of Z2. 



{g G'l^{x,XQ)Gl^{y,XQ)) one gets 
^ 1 



(G^) 1 



(31) 
(32) 



xF-,{x)Al{x)F-,{y)Al{y)D{ax - Py) + 0{D) 

1 

•^2 = .r2^ . / d^^d^y r''F^,{y)A>'Jy)FUx)y, f adaD{x - ay) + 0(^0) (33) 





Thus we have obtained effective action for instanton in NP vacuum in bilocal approximation"^. It can be seen from 
numerical calculation that typical instanton size in QCD is pc 0.25 fm. For such p instanton field in it's center is 
strong F^^{x — xo,pc) — 192/ p* ^ (G^), and therefore classical instanton solution is not much deformed in the area 
I a: I < p, which gives the main contribution to integrals (^-(p3|). 

The problem of asymptotic behavior of instanton solution far from the center |a;| 3> p was studied in detail in 
Refs. H, ||. Our numerical analysis shows that the value of pc is almost not affected by the asymptotic of 
classical instanton solution provided that condensate (G^) and correlation length Tg have reasonable values. 

Next we have to define function D{z). We use Gaussian form D{z) = exp{— fi^ z^) ,u, = 1/Tg. D{z) is a monotonously 
decreasing function with typical correlation length Tg, and numerical calculations ||2q| show that explicit form of D{z) 
is not important. For instance, substitution of exp(— /i^z^) with exp(— /iz) almost does not affect the value of pc- 

We use the standard profile function / = p^ /{p^ +a;^). Certainly, exact instanton profile (i.e. profile that minimizes 
action) is different from this one. Nevertheless, knowing the dependence of ScH on p we will be able to make a 
conclusion about instanton behavior (inflation or stabilization) and to determine it's typical size in NP vacuum. 

Thus we obtain after integration over space directions 

oo 

2^^(G^) ,4 [ x^ 

^'■- = ^1I^nT^I^ 7 (.2 + ^2)2 ^(^) (34) 





^para - ^4 ^4 ^2 _ ^ ^ j "^^"^^ (2,2 + ^2)2(^2 + ^2)2 ^^^i 



00 11 

S. = ^^]^C« / dxdyj-^^^^^^^^ J da J dPe'^-'^'+^'y'\h{2aPxy) + h{2aPxy)) (36) 





where In{x) = e~ ^ Jn{ix) is Bessel function of imaginary argument (Infeld function), = pp and 

1 1 

,(.) =jadaj me-^-^^'^' ^ e-^^ (3^ ' 6^) + 6^ - 2^ + Tx"^^^ ^''^ 


X 2 

where $(x) = J is error function. 





The tensor structure of instanton was used to derive 



We studied asymptotic behavior of instanton effective action and it's dependence on dimensionless parameter 
C = ^p. For a small-size instanton p ^ (C — > 0) we found for 

^dia - -^^^^(G2)/ln(M, (39) 

+ ^1 + ^2 - + O {{G^)p\^^p)) . (40) 



'S'para 



g4 7V2 

In the opposite limit of large p » (C — > oo) one has 



S'para + 5*1 + ^2 ^ COUSt . (42) 



Differential instanton density is proportional to ^j^^ oc exp{—Scf[)- Thus, " diamagnetic" interaction of instanton 
with NP fields leads to its infrared stabilization in size p. 

IV. DISCUSSION AND CONCLUSION 

In this work we have studied instanton behavior in NP vacuum, which was parametrized with gauge invariant 
vacuum averages of gluon field strength. We have derived effective instanton action in bilocal approximation and 
have shown that "diamagnetic" term S'dia leads to IR stabilization of instanton. Numerical results for Scs arc plotted 
in Fig. 1^. Three curves correspond to SU{2) and SU{3) gluodynamics and SU{3) QCD with quarks. Instanton 
size distribution dn/d^zdp ~ exp(— S'eff) and corresponding lattice data |^ are presented in Fig. ^. All graphs are 
normalized to the commonly accepted instanton density 1 fm^*. It should be mentioned that over the last years 
the lattice results have not much converged (compare the recent conference [p8[). Different groups roughly agree 
on instantons size within a factor of two, e.g. p = 0.3... 0.6 fm for S'J7(3) gluodynamics. There is no agreement 
at all concerning the density N/V. As a tendency, lattice studies give higher density and larger instantons than 
phenomenologically assumed. 

Using our model we find for (G^) — 5.92 GeV^ and Tg — 0.22 fm that pc — 0.15 fm, which is less than phenomcno- 
logical 0.3 fm) and lattice results. However, we can present physical arguments to explain these deviations. Lattice 
calculations include cooling procedure, during which some lattice configurations of gluon field are discarded. This 
procedure can result in a change in gluon condensate (G^), and thus instanton size distribution is calculated at a value 
of gluon condensate {G'^)cooi which is smaller than physical value (G^) Therefore, lattice data for average instanton 
size j5 should be compared with our calculations for pc, performed at smaller values of (G^). We show dependence 
of Pc on (G^) for several values of Tg in Fig. ||. One can see that increase of (G^) results in decrease of instanton 
size, and that effect is a result of nonlocal "diamagnetic" interaction of instanton with NP fields. Fig. |^ shows pc as 
a function of Tg for several values of (G^). The smaller Tg, the larger instanton is. It is physically clear that less 
correlated NP fields {Tg 0) have smaller influence on instanton field configuration, which occupies 4D euclidian 
volume with characteristic size p {p ^ Tg). On the other hand, perturbative quantum fluctuations tend to inflate 
instanton, and therefore pc increases with the decrease of Tg. In full QCD (see Table |) we find pc ~ 0.25 fm. 

We did not go beyond bilocal approximation in this work. As mentioned above, this approximation is good 
enough not only for qualitative, but also for quantitative (with accuracy of several percent [^|) description of some 
phenomena in nonperturbative QCD. In the problem under consideration there are two small parameters. These are 
l/g^{pc) ~ 0.15 . . .0.25 and 1/Nc, and there powers grow in each term of cluster expansion. Moreover, we made an 



estimate for the sum of leading terms in cluster expansion 1 26 , and found that IR stabilization stays intact (pc appears 



to be a little smaller). Thus, proposed model describes physics of single instanton stabilization in NP vacuum, not 
only qualitatively, but also quantitatively with rather good accuracy. 



Besides, cooling probably affects correlation length Tg as well. 
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FIG. 2: Instanton density dn/d^zdp and lattice data [ p7{ . 

SU{?i) full (solid line), SU{2) quenched (dotted line) and SU{3) quenched (dashed line) 
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FIG. 3: Instantou size as a function of gluon condensate (JVc = 3, Nf = 4) at Tg = 0.2 fm (dotted line), Tg = 0.3 fm (dashed 
line), Tg = 0.34 fm (solid line) 
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FIG. 4: Instanton size as a function of correlation length {Nc = 3, Nf = 4) at (G^) = 0.5 GeV* (dotted hne), (G^) = 0.87 GeV 
(solid line), (G^) = 1.0 GeV (dashed line) 



